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Abstract 

A  tomographic  array  is  placed  in  a  2-layer,  fiat  bottom,  steady-wind  driven  quasi-geostrophic 
circulation  model  to  investigate  whether  the  analysis  of  acoustic  travel  time  changes  can 
detect  large-scale  barotropic  oscillations.  Time  series  of  sea  surface  elevation  and  upper 
-  and  lower  layer  meridional  currents  are  generated  for  comparison  against  a  series  of  acous¬ 
tic  travel  times.  The  spectra  of  these  time  series  exhibit  a  broad  mesoscale  peak  near  a 
period  of  40  days.  The  spectrum  of  the  acoustic  travel  time  contains  a  significant  peak 
due  to  a  resonant  barotropic  oscillation  with  a  period  of  28.6  days  which  is  not  present  in 
the  spectra  of  the  point  measurements.  In  this  numerical  model,  basin-scale  tomographic 
measurements  are  a  better  method  of  sensing  the  large-scale  resonant  barotropic  oscillations 
than  are  conventional  point  measurements  because  the  tomographic  system  attenuates  the 
“noise”  from  the  mesoscale. 
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Chapter  1 


1.1  Introduction 

The  resonant  barotropic  oscillations  (modes)  predicted  by  the  quasi-geostrophic  (QG) 
theory  are  difficult  to  observe  in  the  ocean.  Large-scale  measurements  over  lengthy  dura¬ 
tions  may  be  required  in  order  to  observe  these  waves. ^or  example,  Luther  (1982)  presents 
evidence  from  2  to  7  years  of  sea  level  records  from  Pacific  islands  that  suggests  the  presence 
of  a  4  to  6  day  barotropic  planetary  oscillation. 

The  capability  of  receiving  long-range  sound  transmissions  in  the  ocean  implies  that 

V  • 

acoustic  methods  may  be  well-suited  for  the  study  of  large-scale  motions.  In  particular, 
basin-scale  tomographic  arrays  may  provide  the  coverage  over  sufficient  time  periods  as 
well  as  the  sensitivity  needed  to  resolve  the  modes  (Bushong,  1987  and  Spiesberger  et  al., 
1988a). 

In  this  study,  we  investigate  whether  the  analysis  of  the  travel  time  fluctuations  of 
tomographic  transmissions  can  detect  large-scale  barotropic  planetary  oscillations  in  a  QG 
modely  We  use  four  years  of  data  from  the  2-layer,  steady-wind,  flat  bottom  eddy-resolving 
general  circulation  model  in  which  Miller  (1986)  detected  b:\rotropic  modes.  A  tomographic 
array  is  placed  in  the  model  and  a  time  series  of  travel  time  changes  is  computed.  We 
compare  this  series  against  series  of  sea  surface  elevation  and  current  to  determine  if  acoustic 
tomography  can  provide  an  effective  method  of  detecting  barotropic  modes. 
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Section  1.2  describes  how  tomographic  travel  times  are  influenced  by  fluctuations  of  the 
temperature,  pressure  and  currents  in  a  QG  model.  Section  1.3  describes  the  response  of 
the  acoustic  travel  times  due  to  resonant  barotropic  oscillations.  Section  1.4  describes  the 
geometry  and  location  of  the  tomographic  section  in  the  QG  model.  Section  1.5  presents 
results  of  the  analysis  of  the  tomographic,  sea  surface  elevation,  and  current  time  series.  A 
discussion  follows  in  section  1.6. 


1.2  The  Forward  Problem 

A.  Definitions 

The  change  in  the  travel  time  along  the  acoustic  ray  path  To  at  geophysical  time  t  is 
approximately  given  by, 

where  the  reference  sound  speed  profile  is  co(*)  and  the  increment  of  ray  path  is  ds  with 
unit  vector  S.  A  Cartesian  coordinate  system  is  adapted  where  x  is  positive  to  the  east  with 
unit  vector  t,  y  is  positive  to  the  north  with  unit  vector  j  and  z  is  positive  upwards  with 
unit  vector  k.  Vectors  are  denoted  by  arrows. 

The  travel  time  is  affected  by  three  processes.  One,  the  vertical  displacement  of  water 
(positive  up),  ij(s*,t),  approximately  changes  the  in-situ  speed  of  sound  by, 

6c{?,t)  g  ri{8,t)-et& ,  (1.2) 

where  the  potential  sound  speed  is  cp  (Munk  et  al.,  1981).  Two,  changes  in  the  hydrostatic 
pressure,  6p(s,t),  which  arise  from  surface  displacements,  alter  the  in-situ  sound  speed  by, 

Sc{a,t)  2  >  (1-3) 

where  the  proportionality  constant,  7,  is  approximately  1.682  x  10~6  m  s-1  Pa-1 
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^0 


z  =0 

z  =  -Ht 


z=-H 


Fig.  1.  Definition  figure  of  the  quantities  computed  in  the  quasi-geostrophic  model. 
The  undisturbed  depth  is  H  and  the  undisturbed  layer  thicknesses  are  Hi  (1000 
m)  and  Hz  (4000  m).  The  interface  displacement  is  h  and  the  surface  displacement 
is  rjo  (both  positive  up).  The  streamfunction  in  the  upper  and  lower  layers  is  ipi 
and  fa,  respectively.  The  densities  in  the  layers  are  pi  and  p3. 


(Munk,  1974).  Three,  the  projection  of  the  horizontal  current, 

u(s,t)  =  u(s,t)i+v(s,t)j  ,  (1.4) 

onto  the  unit  vector,  a,  along  the  ray  path  also  affects  the  travel  time.  The  position  of  a 
ray  is  given  by, 

8  —  Xr ayi  ■+•  J/roy 3  "k  zroy^ 

Vertical  currents  are  negligible  in  this  study. 

B.  Relation  to  the  Quasi-Geostrophic  Model 

The  quantities  computed  in  the  quasi-geostrophic  model  are  shown  in  Figure  1.  The 
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vertical  water  displacement  is  given  by, 

—■teiVti*  z  >  -H1  (upper  layer) 
^  z  <  ~Hi  (lower  layer) 


f){x  ,y,z,t)  =  { 
The  perturbation  pressure  is, 


(1.5) 


Sp{x,y,z,t)  = 


Pofofo{x,y,t) 


z  >  -H\  (upper  layer) 


(1.6) 


Po9Vo{x,y,t)+  pog'h(x,  y,t)  z  < -Hi  (lower  layer) 

where  the  mean  density  is  po,  the  Coriolis  parameter  is  fo  and  the  reduced  gravity  is, 

,  _  s{Pz  ~  Pi) 

9  Po 


The  interface  displacement  is  given  by, 


r/0> 


h{x,y,t)  =  (  — )(^i(x,y,f)  -  fo (x,y,t))  , 

(Holland,  1978).  The  horizontal  currents  are  given  by, 

«((*■*.  0  -  (.  =  1,3)  , 


(1.7) 


(1.8) 


Parameter  values  are  po  =  1025.0  kg  m-8,  fo  =  9.3  x  10~8  s-1  and  g'  =-  .02  m  s~2.  The 
origin  of  the  coordinate  system  is  the  lower  left  hand  corner  of  the  model  (Fig.  5). 

The  QG  model  does  not  include  thermodynamics.  Therefore,  the  potential  sound  speed 
gradient  is  taken  from  a  climatological  average  in  the  Northeast  Pacific  (Fig.  2).  The 
gradient  is  small  below  a  depth  of  about  one  kilometer  because  the  ocean  is  nearly  adiabatic 
in  this  region. 
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POTENTIAL  SOUND  SPEED  GRADIENT  (s  ') 

Fig.  2.  The  climatological  potential  sound  speed  gradient  in  the  Northeast  Pacific 
near  30°N,  150°W. 

1.3  The  Forward  Problem  for  Modes 

Miller  (1986)  demonstrates  that  the  barotropic  oscillations  in  the  QG  model  consist  of  a 
summation  of  resonant  Rossby  (planetary)  basin  modes  which  are  excited  by  the  mesoscale 
turbulence  of  the  free  jet  region.  Solutions  for  the  flat  bottom,  normal  modes  are  given  by, 

V'nm(*,y,0=:R«{-Anmexp(i^-)sin(^)sin(~^)exp(»(rnmt)}  ,  (1.9) 

where  Anm  is  the  amplitude  coefficient  (m*  s-1),  is  the  gradient  of  the  Coriolis  parameter, 
L  is  the  length  of  a  side  in  the  square  basin,  n,m  are  positive  integers  and  Re{x}  denotes 
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the  real  part  of  x  (Pedlosky,  1979).  The  modal  frequency  is  defined  as, 


0 

a"m  ~  ' 

We  seek  to  determine  the  travel  time  changes  due  to  the  individual  modes.  The  travel 
time  changes  due  to  the  pressure  and  temperature  fluctuations  associated  with  these  modes 
are  about  an  order  of  magnitude  smaller  than  the  fluctuations  due  to  currents  (Spiesberger 
et  al.,  1988a).  Equation  (1.1)  simplifies  to, 

6T(t)S£-\f  u-8  ds  ,  (1.10) 

C0  JTo 

where  u  is  the  current  associated  with  the  mode,  and  the  small  variation  of  the  sound  speed 
is  neglected.  The  horizontal  currents  are, 


.  .  .  mx.  /  0X  \  , mjr y.  .  ,  . 

«(*.» »*)  =  ^nm (  j  )  cos( ~~ — )cos(— -)sm(-— )cos(onmt)  , 

(i.ii) 

r  4\  A  frn*\  r  Px  \  tnnx\ 

u(x,y,t)  =  A„^(y)cos(^— )cos(— )  “ 

{^Z — )  sm(~— )sxn(— ))  sin(  -  )  cos(nnmt)  . 

“nm  “nm  L>  ’  L, 

For  a  tomographic  section  oriented  along  a  north-south  direction,  equations  (1.10)  and 
(1.11)  yield, 

fnt)s.^A  ,  (1.12) 

c0 


where  the  range  between  the  tomographic  instruments  is  R. 


1.4  Selection  of  Section  Location 

We  want  to  place  the  tomographic  array  in  a  portion  of  the  QG  model  that  is  repre¬ 
sentative  of  the  Northeast  Pacific  Ocean  because  this  study  is  investigating  the  feasibility 


of  utilizing  e  3000  by  4000  km  tomographic  array  in  the  Northeast  Pacific  to  search  for 
such  modes  (Spiesberger  et  al.,  1988a).  Using  Roden’s  section  of  the  Pacific  along  158°W 
as  a  reference,  typical  eddy-like  features  (wavelengths  of  100  to  200  km)  perturb  the  in-situ 
sound  speed  by  about  0.2  to  0.5  m  s-1  at  1000  m  depth.  This  corresponds  to  vertical 
displacements  between  20  to  50  m  (Fig.  3). 

A  significant  problem  of  this  QG  model  is  the  inadequate  zonal  penetration  of  both 
the  eastward,  mid-latitude  jet  and  of  the  abyssal  mean  and  eddy  kinetic  energies  (Schmitz, 
1982).  The  model’s  deep  gyre  extends  to  1000-1200  km  from  the  western  boundary,  whereas 
Worthington  (1976)  found  that  it  extends  to  2500  km  in  the  North  Atlantic.  As  a  result, 
direct  comparisons  cannot  be  made  between  the  model  and  the  ocean  at  equivalent  distances 
but  a  scaling  factor  must  be  used  instead. 

The  QG  sea  surface  elevation  spectra  levels  are  comparable  to  those  of  Pacific  island 
stations  only  in  the  band  between  25-50  days  period  (Miller,  1986).  On  either  side  of 
this  band,  the  model  spectra  levels  computed  at  a  point  just  beyond  the  farthest  eastward 
penetration  of  the  free  jet  are  roughly  10  dB  lower  than  those  computed  by  Luther  (1982) 
for  several  stations. 

The  resonant  modes  within  this  band  usually  have  longer  wavelengths  in  the  north- 
south  direction  than  in  the  east-west  direction  (see  Table  1;  Miller,  1986).  To  maximize  the 
response  of  the  acoustic  travel  times,  a  north-south  section  iB  desirable.  We  focus  our  efforts 
on  detecting  mode  17,  and  place  a  4000  km  length  section  at  x=1160  km  where  the  amplitude 
of  the  modal  current  in  the  meridional  direction  has  a  maximum  value.  Parameters  for  this 
mode  are  listed  in  Table  1.  At  this  value  of  x  in  the  QG  model,  the  interface  displacements 
(at  1000  m  depth)  are  about  5  to  50  m  which  are  comparable  to  Roden’s  observations 
(Fig.  4). 

The  position  of  the  tomographic  section  with  respect  to  the  mid-latitude  jet  is  shown  in 
Figure  5.  Measurements  of  the  sea  surface  elevation  and  the  meridional  current  are  taken 
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Fig.  3.  A  meridional  section  of  sound  speed  in  the  Northeast  Pacific.  Contour 
intervals  are  at  1  ms-1.  We  thank  Roden  for  permission  to  use  his  figure. 

at  y=1500  km,  and  we  refer  to  this  location  as  station  1.  The  free  jet  meanders  with  time 
and  its  effects  are  felt  approximately  600  km  to  either  side  of  the  mid-basin  point  (y=2000 
km).  On  the  average,  the  section  is  just  within  the  farthest  penetration  of  the  jet,  the  region 
at  which  the  eddies  break  off  and  travel  toward  the  western  boundary.  Station  1  and  the 
central  portion  of  the  tomographic  section  are  influenced  by  eddy  processes. 

In  summary,  it  is  difficult  to  locate  a  region  of  this  model  whose  eddy  field  is  similar 
to  that  found  in  the  Northeast  Pacific.  However,  we  believe  that  the  selected  section  is 
realistic  enough  so  that  a  comparison  of  the  tomographic  and  point  measurements  provides 
a  useful  guide  for  further  investigation. 


(n,  m) 

Period  (days) 

A„  (km) 

A„  (km) 

rms  Amp.  (cm) 

rms  Vel.  (cm  s  x) 

5,1 

29.1 

1600 

8000 

1.1 

.  0.7 

Table  1.  The  scales  of  mode  17  in  this  QG  model.  The  x  and  y  wavelengths  are 
indicated  along  with  the  rms  values  of  the  maximum  sea  surface  displacement  and 
the  maximum  meridional  current. 


1.5  Spectral  Analysis 

Four  years  of  statistically  steady  state  streamfunction  values  are  sampled  at  4-day  in¬ 
tervals.  Three  acoustic  travel  time  series  are  generated  from  three  acoustic  rays  traveling 
between  a  source  and  a  receiver.  The  source  and  receiver  are  at  depths  of  800  m  and  are 
separated  by  4000  km  (Fig.  5).  The  angles  of  the  rays  at  the  source  (with  respect  to  the 
horizontal)  are  11.7°,  6.3°  and  2.1°,  and  their  upper  turning  depths  are  154  m,  447  m  and 
638  m,  respectively.  The  rays  are  traced  through  a  climatological  sound  speed  profile  in 
the  Northeast  Pacific  at  the  same  location  at  which  the  potential  sound  speed  gradient  is 
computed  (Fig.  2).  The  travel  time  series  for  the  11.7°  ray  is  detrended  by  a  smoothed 
cubic  spline  so  that  periods  greater  than  about  one  year  are  removed  (Fig.  6).  Time  series 
are  recorded  for  sea  surface  elevation  along  with  upper  and  lower  layer  meridional  currents 
for  comparison  against  the  series  of  detrended  travel  time  change  for  the  11.7°  ray  (Fig.  7). 
The  fluctuations  of  sea  level  and  meridional  current  have  a  predominant  period  of  40-45 
days.  A  dominant  period  is  not  visually  evident  in  the  travel  time  fluctuations. 

The  power  spectral  density  levels  are  computed  by  dividing  each  time  series  into  4  non¬ 
overlapping  segments,  each  of  which  has  a  duration  of  approximately  one  year.  A  Bartlett 
window  is  applied  to  the  data,  so  the  Bpectral  estimates  are  Chi-squared  random  variables 
with  12  degrees  of  freedom  (Jenkins  and  Watts,  1968).  The  highest  sidelobe  level  for  this 
window  is  -25  dB  (Oppenheim  and  Schafer,  1975). 

Within  the  band  of  interest  (25-50  days) ,  the  sea  level  and  current  spectra  are  similar 


15 


HEIGHT  (m) 


0  1000  MOO  3000  4000 

EAST -WEST 


Fig.  5.  Upper:  Map  of  the  upper  layer  streamfunction  on  day  4.  The  solid  line 
indicates  the  position  of  the  tomographic  section  (x=1160  km)  and  the  asterisk 
marks  the  location  of  station  1  where  the  sea  surface  elevation  and  currents  are 
computed  (y=1500  km).  Contour  interval  is  1  X  104  m*  s-1  and  the  dashed  lines 
correspond  to  negative  values.  Bottom:  Same  except  for  the  lower  layer.  The 
contour  interval  is  0.25  x  104  mJ  s_l. 
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in  shape  (Fig.  8).  The  three  spectra  exhibit  a  broad  peak  centered  near  .025  cpd  (40  day 
period),  and  are  similar  to  Miller’s  results  (see  his  Figures  3b,  4e  and  5e).  The  spectrum  of 
the  travel  time  change  for  the  11.7°  ray  also  displays  a  broadband  peak  in  the  same  band 
along  with  a  peak  at  a  period  of  28.6  days  which  is  significant  at  the  95%  confidence  level 
(Fig.  8).  This  peak  is  also  significant  at  the  95%  level  in  the  spectra  of  the  other  two  rays 
(not  shown). 

Is  this  peak  caused  by  mode  17?  We  tentatively  identify  the  peak  as  mode  17  which  has 
a  theoretical  period  of  29.1  days.  Using  equation  (1.12),  we  calculate  the  travel  time  change 
due  to  this  mode.  If  we  assume  the  modal  energy  is  spread  uniformly  in  a  bandwidth  A / 
(cpd),  the  power  spectral  density  level  P0  (s*  cpd-1)  is  given  by, 


where  cr|f  is  the  variance  of  the  travel  time  fluctuation  caused  by  the  mode.  Table  2  lists 
the  predicted  and  measured  rms  travel  time  change  and  spectral  level  for  this  mode.  The 
predicted  spectral  level  of  3.18  x  10-i  (ss  cpd-1)  is  computed  using  a  half-power  bandwidth 
of  .002  cpd  which  is  a  Q  of  about  17  (Miller,  1986).  The  observed  spectral  level  of  this  peak 
is  about  1.69  x  10-J  (s*  cpd-1).  The  predicted  and  observed  levels  are  similar. 

We  seek  to  determine  the  influence  of  the  eddy  field  on  the  travel  times  to  ensure  that 
the  peak  is  not  the  result  of  mesoscale  processes.  The  QG  model  is  high-passed  filtered  using 
a  Bartlett  window.  For  this  filter,  wavelengths  of  800,  1000  and  2000  km  are  attenuated 
by  -3.4,  -6.1  and  -16.0  dB,  respectively.  Therefore,  the  mesoscale  field  is  passed  and  the 
long-scale  planetary  modes  are  attenuated.  The  averaged  periodogram  of  the  travel  times 
is  computed  over  the  four-year  record  for  the  high-passed  model  (Fig.  9).  The  spectral  level 
is  about  30  dB  less  than  the  peak  at  28.6  days  observed  in  the  unfiltered  model.  Therefore 
we  conclude  that  the  peak  at  a  period  of  28.6  days  is  caused  by  a  planetary  basin  mode 
rather  than  the  mesoscale  field. 
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Fig.  7.  Time  series  of  sea  surface  elevation,  upper  layer  meridional  current,  lower 
layer  meridional  current  and  detrended  acoustic  travel  time  change  for  the  11.7° 
ray  (from  top  to  bottom).  Elevation  and  current  are  measured  at  x=1160  km, 
y=1500  km  (Fig.  5).  The  first  2  years  of  the  4-year  records  are  displayed.  The 
values  are  computed  at  4-day  intervals. 
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rms  Travel  Time  Change  (ms) 

Spectral  Level  (s2  cpd  *) 

Predicted 

-11.27 

3.18  x  10"2 

Measured 

7.82 

1.69  x  10~5  “1 

Table  2.  The  predicted  values  of  the  rms  travel  time  change  and  the  power  spec¬ 
tral  density  level  are  compared  with  the  values  computed  from  the  tomographic 
array.  The  predicted  rms  travel  time  change  is  obtained  from  eqn.  (1.12). 


1.6  Discussion 

Although  basin-scale  resonant  barotropic  modes  have  been  detected  in  numerical  models 
(for  example  Willebrand  et  al.,  1980;  Miller,  1986),  attempts  to  detect  this  type  of  motion 
in  the  real  ocean  have  not  been  as  successful.  Luther  (1982)  presented  the  most  convincing 
but,  by  his  own  admission,  still  inconclusive  evidence  of  a  barotropic  vorticity  mode  from 
his  analysis  of  sea  level  records  covering  a  wide  geographic  area  in  the  Pacific. 

Luther  (1983)  gives  two  reasons  for  the  lack  of  direct  observations  of  basin-scale  vorticity 
modes  in  the  ocean.  First,  because  the  modes  are  weakly  excited,  their  signal  may  not  be 
strong  enough  relative  to  the  background  noise  to  be  extracted  from  the  data  he  used. 
Second,  topography  reduces  the  period  and  scale  of  a  planetary  mode,  changing  it  into  a 
topographic  mode  that  is  restricted  to  the  vicinity  of  the  controlling  feature. 

In  this  QG  model,  two  factors  hinder  the  observations  of  barotropic  planetary  modes  in 
the  sea  surface  and  current  records.  First,  the  mesoscale  field  masks  the  signals  of  the  modes 
in  the  point  measurements.  Second,  because  the  modes  are  closely  spaced  in  frequency,  the 
presence  of  friction  and  non-linearities  can  cause  the  narrow  modal  peaks  to  broaden  into  a 
broadband  peak  (Miller,  1986).  However,  the  tomographic  system  is  not  as  strongly  affected 
by  these  factors.  A  basin-scale  tomographic  array  attenuates  the  signals  from  the  mesoscale 
processes  so  that  a  mode  near  28.6  days  can  be  detected  at  the  95%  confidence  limit  for  a 
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Fig.  8.  Power  spectral  densities  of  sea  surface  displacement  (upper  left)  and 
meridional  currents  in  the  upper  and  lower  layers  (lower  left  and  upper  right)  at 
station  1  and  acoustic  travel  time  change  for  the  11.7°  ray  (lower  right)  computed 
over  the  4-year  record.  In  the  travel  time  figure,  the  period  and  rms  travel  time 
of  mode  17  are  labeled.  In  the  other  figures,  the  arrow  marks  the  period  of  the 
mode.  The  95%  confidence  limits  are  shown. 
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Fig.  9.  The  power  spectral  density  of  acoustic  travel  time  change  for  the  11.7®  ray 
traveling  through  the  eddy  held  (but  not  the  large-scale  field).  The  arrow  marks 
the  period  of  mode  17,  and  the  95%  confidence  limit  is  shown. 
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record  length  of  four  years.  In  other  words,  the  acoustic  travel  time  changes  are  given  by  a 
spatial  integral  between  the  source  and  the  receiver  which  preferentially  attenuates  the  short 
wavelengths  in  the  ocean  (Spiesberger  et  al.,  1988a).  Also,  the  tomographic  measurements 
are  able  to  resolve  an  individual  mode  by  placing  the  section  at  a  position  where  the  modal 
current  has  a  maximum  value. 

Spiesberger  et  al.,  (1988a, b)  have  demonstrated  that  the  acoustic  phase  of  the  tomo¬ 
graphic  signals  might  be  used  to  make  estimates  of  relatively  slow  barotropic  currents 
whose  length  scales  are  large  (order  of  1000  km  or  more).  No  one  knows  what  the  spectra 
of  acoustic  travel  times  over  basin-scale  distances  look  like  for  periods  greater  than  a  few 
days.  However,  the  internal  wave  field  imposes  a  limit.  For  a  tomographic  array  consisting 
of  one  source  and  ten  receivers  at  about  4000  km  range,  the  threshold  at  which  a  basin- 
scale  current  can  be  detected  with  a  confidence  of  95%  is  given  by  uS  33.6T  micro-meters 
per  second  where  the  period  of  the  oscillation,  T,  is  expressed  in  days.  For  mode  17,  T  is 
approximately  29  days,  and  the  threshold  speed  is  0.1  cms"1.  The  amplitude  of  the  mode 
17  current  in  this  QG  model  is  1.0  cm  s_1.  Therefore,  if  internal  waves  limit  the  precision  of 
the  tomographic  array  and  the  real  ocean  contains  modes  whose  scales  are  similar  to  those 
in  the  QG  model,  then  the  modes  would  be  detected  by  the  tomographic  system. 

Platzman  (1981)  has  computed  the  periods  of  the  ocean’s  modes  for  the  range  between  8 
and  80  hours.  Any  experimental  program  which  searches  for  modes  with  periods  exceeding 
80  hours  would  be  greatly  enhanced  by  an  extension  of  the  numerical  calculations  to  longer 
periods.  However,  at  these  longer  periods,  the  stratification  of  the  ocean  and  the  non¬ 
linear  interaction  of  the  modes  with  other  large-scale  currents  near  the  same  period  may  be 
important  both  for  the  interpretation  of  the  data  and  for  a  reliable  numerical  computation 
of  the  modal  structure. 
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